Electron-radiation interaction in a Penning trap: 
beyond the dipole approximation 
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We investigate the physics of a single trapped electron interacting with a radiation field without 
the dipole approximation. This gives new physical insights in the so-called geonium theory. 

PACS number(s): 03.65.Bz, 42.50.Dv, 12.20.-m 



I. INTRODUCTION 



Simple systems like single electron or ion provide a useful tool to investigate the fundamental laws of nature. Hence 
in the past decades there has been increasing interest on trapping phenomena |jj . It is now routinely possible to trap 
a single ion [Q, and it would allow us to study QED (Quantum Electro- Dynamics) also when the trapped ion interacts 
with a radiation mode. On the other hand the electron stored in a Penning trap |j| permitted accurate measurements 
\ too Q . This system has been called " geonium atom" since it resembles an hydrogen atom where indeed the binding 
for the electron is to an external apparatus residing on the earth ||). The geonium system was recently studied to 
implement some interesting quantum optics situations, like QND (Quantum Non-Demolition) measurements Q and 
the generation of nonclassical states 0. 

Here, just after one hundred years from the electron discovery, we would show new interesting quantum features 
. arising on a trapped electron interacting with the radiation field, when no dipole approximation is made. 
ON \ It is well known that in the geonium system Q the motion of the electron can be separated into three independent 
pH , harmonic motions: the axial, cyclotron and magnetron. On the other hand, it is also well established that entangled 
O systems are extremely interesting for many purposes. 

In the present work we propose a way of coupling the three harmonic oscillators of the geonium system by simply 
superposing a radiation field to the trapping fields. More concretely, we show that when the trapped electron os- 
cillates in a standing wave field, there could be linear, or nonlinear coupling among the axial motion and the other 
ones; although, in particular we will only consider the axial-cyclotron interaction. Hence, we shall present the more 
immediate consequenses of such an entanglement, like indirect measurements on the cyclotron mode, then we shall 
investigate the generation of nonclassical features. Moreover, the analysis in all cases will be performed by taking into 
account the environmental effects as well. 

The paper is organized as follow: Section II is devoted to the description of the model. The first order coupling 
(linear) between axial and cyclotron motion is considered in Section III, while in Section IV the second order coupling 
is discussed. In Section V we further discuss the possibility of generating nonclassical states. Finally, we present 
conclusions. 



II. THE MODEL 

The geonium system consists |j| of an electron of charge e and mass too moving in a uniform magnetic field B, 
along the positive z axis, and a static quadrupole potential 

where d characterizes the dimension of the trap and Vo is the potential applied to the trap electrodes JbJ . 

In this work, in addition to the usual trapping fields, we embed the trapped electron in a radiation field of vector 
potential A ex t . To simplify our presentation, we assume the a priori knowledge of the electron spin [^) . Neglecting 
all spin-related terms, the Hamiltonian for the trapped electron can then be written as the quantum counterpart of 
the classical Hamiltonian 
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where c is the speed of light, and 



A = -f A B + A ext . 



(3) 



where f = {x,y, z) and p = (p x ,p y ,p z ) are respectively the position and the conjugate momentum operators of the 
electron. 

The motion of the electron in absence of the external field A ex t is the result of the motion of three harmonic 
oscillators ||: the cyclotron, the axial and the magnetron, which are well separated in the energy scale (GHz, MHz 
and kHz respectively). This can be easily understood by introducing the ladder operators 
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where the indexes z, c and m stand for axial, cyclotron and magnetron respectively. The above operators obey the 
commutation relation [cti,ctj] = 1, i = z, c, m. The angular frequencies are given by 



eV _ eB 

So, when A ext = 0, the Hamiltonian (Eh simply reduces to 



2uj c . 



1 



H = huj z ( a\a z + ^ I + TiLj c ( a\a c + ^ ) - Tiuj m ( al n a m + ^ 



1 



(7) 



(8) 



Instead, when the external radiation field is a standing wave along the z direction (with frequency f2 and wave 
vector k) and circularly polarized in the x — y plane [pi, we have 



A ext = (-i [ae l!2f - a*e~' nt ] cos(fcz + 0), [ae lUt + a*e" int ] cos(kz + <^),0) 



(9) 



In such a case, and for frequencies 11 close to uj c , we can neglect the slow magnetron motion, and the Hamiltonian 
a) becomes 



H = hu! z ala 
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and the phase (p is the phase of the applied radiation field (i.e. arga). The other phase </> defines the position of the 
center of the axial motion with respect to the standing wave. The third and fourth terms in the right hand side of the 
Hamiltonian ( JToj ) describe the nonlinear interaction between the trapped electron and the standing wave which gives 
rise to a coupling between the axial and the cyclotron motion, whose effect will be analyzed in the next sections. In 
the usual Penning traps the quantity k(z) can reach values up to w 0.1 [||, when Q w lu c . This leads us to explore 
the physics beyond the usual dipole approximation for the cosine term in Eq. (|l0|). The cosine factor cos(fcz + 4>) can 
be split as 



cos(fcz + cj)) = cos <j) cos(/cz) — sin <f> sin(fci) . 



(12) 
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and two typical situations corresponding to 4> — and cj) — ir/2 can be easily exploited. By making the usual dipole 
approximation these two cases correspond to a mere driving term on the cyclotron motion ((f) = 0) or to no effect at 
all = tt/2). 

In the following Sections the behaviour of the trapped electron in these two paradigmatic limits is studied. All the 
other possible values of cj) will give rise to combinations of these two cases and can be easily studied. 

We further note that the last term in ( JlC| ) can be neglected since the parameters ( pl| ) are such that x/\ e \ ~ l e l/ w c- 

III. THE CASE OF <f> = tt/2 

In this section we consider the case <f> = n/2. Developing sin(fcf) in power series and keeping only the first order 
term we can approximate the Hamiltonian (|10|) by 

H = Huj z (^a\a z + ^ + {^cPc + ^ 

+ h [e*a c e int + ea£e~ iat ] kz . (13) 
In the case of perfect resonance, Q — uj c , and in a frame rotating at that angular frequency we get the solution 

z(t) = |z(0) - \e\kxj cos(oj z t) + —p z (0) sin(w,t) + \e\kX v (14) 

L J muj z 

Pz(t) = Pz(0) cos(u> z t) ~ muj z z(0) - V2\e\kX v sin(uj z t) (15) 
where we have introduced the cyclotron quadrature 



Equation ( ^5[ ) suggests us an indirect way to determine the probability distribution for the cyclotronic quadrature, 
V(X V ). We recall that in the geonium system the measurements are performed only on the axial degree of freedom due 
to the non existence of good detectors in the microwave regime. The oscillating charged particle induces alternating 
image charges on the electrodes, which in turn cause an oscillating current to flow through an external circuit. The 
current will be proportional to the axial momentum p z , hence a measurement of this current will also give the value 
of the quadrature X v . Measurements when the standing wave is 'off' should be done preventively to set the initial 
conditions. Then, repeated measurements lead to the desired statistics V(X V ). 

If the procedure is further repeated for several values of the phase if, we obtain the set of marginal probabilities 
V(X, if), which allows the tomographic imaging of the quantum state of the cyclotron mode 

We now consider the effects of the thermal damping through the resistance of the external circuit connected with 
the measurement apparatus. In such a case the equations of motion for the axial degree of freedom become 

at mo 

% = -Lu 2 z m z- ^p z -V2hk\e\X v +Z, (18) 
at mo 

where the noise term is that of Johnson noise with expectation values = 0, and (£(£)£(*')) = 2j z KBT5(t—t'), 

the damping constant j z is proportional to the readout resistor, Kb is the Boltzmann constant and T the equilibrium 
temperature. 

By using the Fourier transforms, we immediately obtain 

c<j z — cjj — iw7 z /mo 

hence the correlation 

WMU _ U)) . gMSfeMi^ . m 

Eq. (|2(]) imposes some limits to the observability of nonclassical effects on the cyclotron motion; in fact the added 
thermal noise should be much less than the cyclotron vacuum noise for the chosen frequency, i.e. ^ z KbT « (7ifc|e|) 2 . 
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IV. THE CASE OF <f> = 



Let us now consider the case of (f> = 0, and keeping only terms up to the second order in kz, the Hamiltonian ( |io| ) 
reduces to 
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which clearly shows the nonlinear coupling as a consequence of the higher order expansion with respect to the case of 
the previous Section. 

Let us study the general case including losses. The latter are present in the axial degree of freedom once the 
connection with the external circuit is established, as pointed out in the previous Section. Instead, the noise on the 
cyclotron degree of freedom could arise e.g. from radiative damping (though it can be strongly reduced with an 
appropriate trap geometry). 

Hence, by starting from the Hamiltonian (pl|), we obtain the following Quantum Stochastic Differential Equations 
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where A — lu c — fi, / is a driving term for the axial motion, j c the cyclotron damping constant, and a™, S are the 
noise terms (we shall consider the situation where only the vacuum contributes to the cyclotron noise). We have 
introduced the scaled variables Z — y/motJz/hz, P z = y/l/hrnou>z p z , 2 = y/lfhrnoUJ z £, and k 2 = Tik 2 /2toow z . From 
the Eq. ( pB| ) we can see that the cyclotron quadrature causes a shift of the resonant frequency of the axial motion, so 
its indirect measurement results feasible. 

The system of equations ( ^2|]25| ) can be linearized around the steady state O]. The stationary values a c , Z and 
P z can be obtained from the following equations 
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The linearized system is then 
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where now the operators indicate the quantum fluctuations with respect to the steady state, and 
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The spectral matrix can be calculated as 
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S(w) = (iujl - M) -1 D(-kjI - M T ) _1 , 



(32) 



where I is the four by four identity matrix, M T means the transposed, and 
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with Nth = KbT/UlOz the number of thermal excitations. 

The momentum correlation for the axial motion will be S44; this quantity is plotted in Fig.l. The dashed line 
represents the resonance in absence of coupling and the solid line the resonance in presence of it. The separation 
between peaks is proportional to the cyclotron quadrature amplitude. So, it gives us an indirect value of that cyclotron 
observable. 

Furthermore, the variance for the amplitude cyclotron quadrature is given by integrating the quantity Sn + S22 + 
S12 + S21, and it is plotted in Fig. 2 (dashed line). The same figure also shows the variance for the orthogonal 
quadrature (solid line). It can be seen that the system exhibits squeezing effects depending on the detuning. It is 
worth noting that such effects are not much sensitive to thermal noise. 

The stability of the system, for the values of parameters used, is checked through the signs of the eigenvalues of the 
matrix M. 

In this Section and in the previous one we have shown that the terms beyond the dipole approximation could play 
an important role and should not be neglected abruptely. As a matter of fact we have presented a variety of effects 
(see e.g. Figs.l and 2) that could be measured in common Penning traps. 

To go further, in the following, we shall explore other possibilities. 

V. NONCLASSICAL STATES 

We now demonstrate the generation of nonclassical effects due to the nonlinearity induced by Hamiltonian (|2 

A. The central resonance 

If we tune the standing wave at frequency fl = lu c , and pass to the interaction picture, the Hamiltonian ( pl| ) simply 
becomes 



H = V2h\e\X c 



1 — k 2 [ ata x H — 



(34) 



where we have disregarded the rapidly oscillating terms a\e~ 2luJ * t and a z e 2%UJzt (i.e. we made the Rotating Wave 
Approximation) . 

Starting from initial coherent states for both modes 

|*(0)> = \a) c ® \p) z , (35) 



we obtain from the Hamiltonian ( 34 ) the following state at the time t 
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(*)> = e-W/ a £-7=flMc®|n>,, (36) 



where 6 n = ien 2 n. In writing the state ( pq ) we have disregarded, for the sake of simplicity, the quantity a— ie(l— n 2 /2)t, 
which is common to each cyclotron component (this corresponds to an overall displacement in the cyclotron phase 
space). 

Therefore, the electron motion evolves classically as a mixture of coherent states. Thus, during the evolution, no 
nonclassical states of the electron are generated. However, because of the entanglement between the cyclotron and the 
axial degrees of freedom, it is possible to generate nonclassical states of the cyclotron motion by performing conditional 
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measurements on the axial degree of freedom. In particular, a measurement of the axial current corresponds to the 
projection onto an eigenstate \p z ) of the axial momentum 



\y(t)) after =Mj2 
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where Af is a normalization constant and {p z \n) z are the harmonic oscillator wave functions in the momentum space. 
It is immediately seen from the above expression, that after the measurement the system is left in a superposition of 
coherent cyclotronic states which could have nonclassical features. Whenever the number of coherent states that are 



being superposed is small, the states are known as Schrodinger cat states 12 . 

It is worth noting that the separation between the superposed coherent states is given by |e|K 2 t, and therefore it 
can be made truly macroscopic emphasizing the nonclassicality (by simply requiring that \e\K 2 t > 1). However, one 
has to be careful when satisfying the above condition, since it also implies a strong excitation of the cyclotron motion 
(the overall displacement that has been disregarded), which in turn could give rise to instabilities or even, the loss of 
the particle over the trap's walls. 

The Wigner function of the cyclotron state generated by conditional measurement can be written as 



W(Q, P) = N 2 J~) c* m c n exp 
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where the variable Q, P are associated to the quadratures X v= q and X v=7r /2 respectively, and 



Cn — @n t • 



(39) 
(40) 



In Fig. 3 we present the Wigner function of the cyclotron state generated by conditional measurement on the axial 
degree of freedom. The negative parts and several oscillations show the nonclassicality of such a state. 

We have considered the measurement process conditioning the cyclotron state as instantaneous, however, it always 
takes a finite time during which the system undergoes the back-action of the measurement apparatus. To take into 
account these effects we should adopt a precise Hamiltonian model describing the measurement of the observable p z . 
Nevertheless, from a phenomenological point of view, we can model the measurement process, performed on the state 
p(t) = \^(i))(^(t)\, during a time r, as the transition 



W + t) \ Pz ){p z 



(41) 



where p(t + r) is obtained from p(t) through free evolution in a thermal bath (representing the back-action of 
the measurement apparatus on the system), while the projector indicates the output resulting at the end of the 
measurement |l3). 

To evaluate the effects of the measurement on the cyclotron state we should calculate the reduced density operator 
(r.h.s. of Eq.(f4l|)). Its corresponding Wigner function is derived in Appendix A as 
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with LJ^ the associated Laguerre polynomials, and T = j z /mo the effective axial damping constant. 

The Wigner function ( f42| ) is plotted in Fig. 4 and shows the deleterious effects of finite time measurement on the 
nonclassical state represented in Fig. 3. Of course, these effects strongly depend on the number of thermal excitations 
N t h as well. 

Once the cat states are generated by the conditional measurements, it would be possible to detect them by using 
indirect measurements as proposed in the previous Sections. 



B. The sideband resonance 



We now return to the Hamiltonian ( J2l| ) to consider another resonance, in this case Q = {oj c — 2u z ) — 5, where 5 is 
a small detuning (i.e. 5 << lu z ) introduced for convenience. In a frame rotating at frequency ui c — <5, we then have 



H = Ti8a],a c — h 



(44) 



This is a trilinear Hamiltonian analogous to that studied in nonlinear optical processes like parametric oscillation or 
second harmonic generation 11 1. 
The equations of motion are 



da c i • 1 i i 2-2 

— = -iSa c + i-\e\K a z , 



da z 
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and, by adiabatic elimination of the cyclotron mode, we get 
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This equation corresponds to an effective Hamiltonian for the axial motion of the type 
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(45) 
(46) 

(47) 
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which shows a well known Kerr-typc nonlincarity. Hence, we should expect nonclassical effects, such as Schrodinger 
cat states, when one starts from the initial conditions ([35]), also in the axial mode. In fact, the evolved axial state can 
be written as 



\^{t)) z = exp [iG {{d\d z ) 2 - d\a z ) t] 



G = 



45 



(49) 



It is easy to show that after a time t — tt/ (2G) the initial coherent state evolves into a cat state of the type discussed 
in Ref. Ill 



= n/(2G))) 



1 

V2 



(50) 



That state shows interference in the phase space which could be detected by measuring an appropriate quadra- 
ture. Therefore, by adjusting the initial conditions we may exploit the axial momentum measurement to see such 
interference. The Wigner function of the state (|50|) results 
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(51) 
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and is represented in Fig. 5. The fact that only two coherent states are being superposed is evident from the two hills 
besides the central structure, differentely from the situation of Fig. 3 where more components contributes to the cat 
state. 

Of course, we should deal again with the problem of measurement, whose process renders the system open, hence, 
the dissipation tends to wash out the nonclassical effects. To evaluate this phenomenon we assume to switch off the 
nonlinearity at the time of cat generation, and a subsequent free evolution of the axial degree of freedom in a thermal 
bath, representing the effects of the external readout circuit. If the latter takes a time r, we have (see Appendix B) 
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(58) 
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(60) 



The Wigner function (|52|) is plotted in Fig. 6 and shows that the cat state fl50| ) is very sensitive to the noise induced 
by the measurement. 



VI. CONCLUSIONS 



In conclusion, we have studied a trapped electron interacting with a standing radiation field and have shown that 
several interesting features can arise when the dipole approximation is not invoked. First, the proposed model provide 
a method for indirect measurement on the cyclotron degree of freedom. On the other hand, the possibilities to generate 
nonclassical states could be useful to test the linearity of Quantum Mechanics , as well as to probe the decoherence 
of a mesoscopic system fl7j] . Furthermore, it is worth noting that the entanglement induced by the radiation field 
could also be used to explore the quantum logic possibilities of a trapped electron system. 

Hence, the geonium system in such configuration could result alternative and/or complementary to other trapped 
systems. In addition it has the advantage of involving a structureless particle, while for example an ion in a Paul 
trap behaves as a two-level system only ideally. Moreover, having the electron an antiparticle, the model could also 
be used to perform some fundamental tests of simmetry. 

Finally, based on theese considerations we conclude that it should be an interesting challenge to experimentally 
implement this model. The realistic values of parameters (see e.g. Ref. |]) we have used yield that feasible with the 
actual technology. The main problem could be represented by the low values of N t h in Sec. IV, however, to better 
evidenciate the desired effects one could adjust the experimental set up in order to increase the dishomcgeneity of the 
field experienced by the particle (to this end, we note that traps bigger than the usual are available as well ]l^|). 
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APPENDIX A 



We consider the position space matrix elements of the state pit), i.e. 

c (Q + Y\ z (Z'\p(t)\Z") z \Q-Y) c 



(61) 



and we denote them as p(Z' , Z") since the evolution will take place only in the axial space. The dependence on the 
cyclotron variables remains implicit. Then, with the aid of Eq. ( |36| ) we get 



p(Z', Z") = ]T C m C* n exp [- (Z' 2 + Z" 2 ) /2\ H m (Z 1 ) H n (Z") , 



(62) 



where H m are the Hermite polynomials, and 

C m = 

c: = 



i\ 1/4 i 



,-1/^/2 Z 



ttJ \j2 m m\ 
l\ 1/4 1 



--(Q + Y\( m ), 



-\p\*/2 (P*) n (c l0 _ Y ) 



(63) 
(64) 



The master equation for the free evolution in a thermal bath ||15| has the corresponding partial differential equation 
for the probability p 

\ i f d 2 d 2 \ T ( d d 

d T p(Z , Z , r) = - - —s j - ~(Z Z) 



rN th (Z'-Z") 2 }p(Z\Z",r), 



(65) 



where we have set T = 7 z /mo. Both, the damping constant and the time are scaled by the axial frequency, i.e. 
T/lu z — > r and tlo z — > t. 

The differential equation (pq) is considerably simplified by the change of variables 



Z = u + v , 
Z" = u-v, 



leading to 

d T p(u,v,T) 
By using the Fourier transform 



i d 2 



2 dudv 



Tv—-4TN th v 2 }p(u,v,T). 



p(u,v)= / dqe 2t " u p(q,v) 



Eq. 



becomes 



dp , , . „ ,dp 2 
— + (q + Tv) — = -ATN th v p , 
OT ov 



which can be solved by the method of characteristics. The solution takes the form 



P(q, v, t) = p (q, (v + ^je rr - ~ 



x exp < - 2N t h v 



1-e 



-2Tt 



(66) 
(67) 



(68) 



(69) 



(70) 



(71) 
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In our case, from the Eqs. fl62[), (|66j), (|67]), (|69|), results 
p(g, 0) = Vttcxp [-w 2 - gr 2 ] 



| £ C m Q2« m ! (-„ - zg)— [2 (, 2 + g 2 )] 

J] |C„| 2 2"n!i m [2 ( V 2 +<7 2 )] 

m— n 

£ C m C:2-n! (-« - iq) m - n [2 (t, 2 + g 2 )] 1 , 

rn>n J 



(72) 



where L™ indicates the associated Laguerre polynomials. Therefore, starting from the above expression, the solution 
( fzT| ) can be easily constructed. 

The Wigner function of the cyclotron state after a measurement giving the result p z (or equivalently P z ), will be 



W(Q, P, t) =Af 2 J dY C (Q + Y\ z (P z \p(t + T )\P z ) z \Q -Y) c e 



2iPY 



(73) 



where TV" is the normalization constant needed after the projection. By inserting identities in terms of the set of states 
{\u ± v) z }, and with the aid of the Fourier transform (|69|), we get 



W(Q, P, t) = TV 2 dY dvp(0,v,r) 



(74) 



The dependence on the cyclotron variables Q and Y is implicit on p. Hence, by performing the integration one arrives 
at the expression (Eg). 



APPENDIX B 



If r is the duration of the measurement, at the end of the measurement we have 



W(Z, P„, T ) = - I dv(Z + v\p z ( T )\Z - v)e 

7T ' 



-2iP« 



dv p(Z, v, r)e 



-2iP z v 



where p(q,v,r) is the same of Eq. (|7l|), but with the initial condition determined by Eq.(|5 

p(q,v,0) = ±e-M 2 -<?-» 2+ ^ 2+ F 2 / 2 

x | exp 2Im(/3) 2 - 2%/2Tm(/3) (? - 2V2iRe((3)v 

+ exp [2Im(/3) 2 + 2y/2iJm(/3)q + 2v / 2£Re(/3> 
— i exp -2Re(/3) 2 - 2V2Re(/3)q + 2\/2Im(/?)v 

+ i exp [-2Re(/3) 2 + 2\/2Re(/?)q - 2 V2Im(/3)u 



(75) 
(76) 
(77) 



(78) 



Thus, performing the double integral in Eq.(f77|) we get the expression (|5 
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FIGURE CAPTIONS 

Fig.l Spectrum of axial momentum for A = 1.5 x 10 4 s -1 , n 2 = 10 -6 , j c = 1.5 s _1 , 7 2 /m = 20 s _1 , |e| = 1.4 x 10 4 , 
s _1 , tp = 37r/4, / = 10 11 s _1 , and Nth = 10 3 - The peak on the right represents the resonance in absence of coupling. 
The separation between peaks is proportional to the cyclotron quadrature amplitude. 

Fig. 2 Variance for the cyclotron quadratures X v —o (dashed line) and X v= ^/ 2 (solid line) as function of the detuning 
A. The values of other parameters are the same of Fig.l. 

Fig. 3 The Wigner function of Eq. fl38| ) is plotted for the parameters 0=1, tn 2 t = — 2Ai, after an axial momentum 
measurement yielding the most probable value of p z . 

Fig. 4 The same of Fig. 3 including the effects of finite time measurement. Here Tt = 0.1 and Nth = 10. 

Fig. 5 The Wigner function of cat state J50| ) is plotted for (3 = 2. 

Fig. 6 The same of Fig. 5 including the effects of finite time measurement. Here T = 6, r = 0.4, and Nth = 10. 
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Fig.l A. M. Martins et al . 
Electron-radiation 
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Fig. 3 A. M. Martins et al . 
Electron-radiation 




Fig. 4 A. M. Martins et al . 
Electron- radiation 




Fig. 5 A. M. Martins et al. 
Electron-radiation 




Fig. 6 A. M. Martins et al. 
Electron-radiation 



